Universal Quantum Computation with Continuous- Variable Cluster States 
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We describe a generalization of the cluster-state model of quantum computation to continuous- 
variable systems, along with a proposal for an optical implementation using squeezed-light sources, 
linear optics, and homodyne detection. For universal quantum computation, a nonlinear element is 
required. This can be satisfied by adding to the toolbox any single-mode non-Gaussian measure- 
ment, while the initial cluster state itself remains Gaussian. Homodyne detection alone suffices to 
perform an arbitrary multi-mode Gaussian transformation via the cluster state. We also propose an 
experiment to demonstrate cluster-based error reduction when implementing Gaussian operations. 
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Introduction — One-way quantum computation [l| 
provides the ability to perform universal quantum com- 
putation (QC) using only single-qubit projective mea- 
surements, given a specially prepared and highly entan- 
gled cluster state. This is in contrast to the traditional 
circuit model, where unitary evolution and coherent con- 
trol of individual qubits are required |2( . Apart from its 
conceptual importance, the cluster-state approach can 
also lead to practical advantages. For example, the re- 
sources required for QC using linear optics 3] can be 
significantly reduced by first creating photonic cluster 
states via nondeterministic gates Recently, a 

four-qubit cluster state has been demonstrated optically 
in the single-photon regime [jj. 

While qubits are typically used in QC, Lloyd and 
Braunstein || proposed the use of continuous variables 
for QC and proved that only a finite set of continuous- 
variable (CV) gates are needed for universal QC. In the 
CV approach, the continuous degree of freedom may be 
used directly, or lower-dimensional systems may be en- 
coded within the modes, such as in the GKP proposal 0, 
which encodes one qubit into each mode. This allows, for 
instance, for the application of standard qubit protocols 
to CV systems. The optical modes of the electromagnetic 
field provide an experimental testbed for these ideas 0] . 

In this Letter, we describe a model of universal QC 
using CV cluster states. We also propose an optical im- 
plementation of our scheme where squeezed light sources 
serve as the nodes of the cluster. The main advantage 
of this approach is that not only can computations with 
the cluster be performed deterministically, but also the 
preparation of the cluster state, including connecting the 
nodes, can be done unconditionally. This is in contrast 
to the discrete- variable linear-optics schemes 0, 0, ^3 > 
where cluster states are created probabilistically. There- 
fore, the CV approach appears to be particularly suited 
for further experimental demonstration of the general 
principles of cluster-state QC. 



In our optical implementation, once the cluster state 
has been created, single-mode homodyne detection alone 
will allow for any multi-mode Gaussian transformation 
to be performed on the information contained within the 
cluster. Analogously to the implementation of Clifford 
gates using qubit clusters, the homodyne detections can 
be done in any order, a property known as parallelism. 
For universal QC, in addition, only one single-mode non- 
Gaussian projective measurement (e.g., photon counting) 
is required. However, parallelism no longer applies to 
non-Gaussian measurements, because the choice of sub- 
sequent measurement bases will depend on the outcome 
of earlier measurements. This adaptiveness of the mea- 
surement bases is again analogous to the qubit case when 
computing non-Clifford gates. While CV cluster states 
have been described previously in Il2j , it is claimed there 
that such states are an insufficient resource for universal 
QC because of their Gaussian character 0] . In fact, they 
are sufficient as long as we can perform a non-Gaussian 
measurement. An analogous result holds for qubit clus- 
ter states, which can be created entirely using Clifford 
group operations |l4| but are nevertheless universal once 
a non-Clifford measurement is allowed. 

Although CV cluster states can be built deterministi- 
cally, it will be impossible to create perfect CV cluster 
states due to the finite degree of squeezing obtainable in 
the laboratory. This results in distortions to the quantum 
information as it propagates through the cluster state. 
We discuss these distortions (along with other errors) 
and propose an experiment that demonstrates how par- 
allelism and post-selection can be used to mitigate these 
effects when implementing Gaussian operations. 

Continuous- Variable Cluster States — Other au- 
thors 0, 0, have extended the cluster-state for- 
malism to d-level systems (qudits). Here we generalize 
these results to CVs. Our use of CVs for QC follows 
the standard prescription given in [J^. The Pauli X 
and Z operators are generalized to the Weyl-Heisenberg 
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group, which is the group of phase-space displacements. 
For CVs, this group is a Lie group with generators q 
and p. These operators satisfy the canonical commuta- 
tion relation [q,p] = i (with K = 1) and when exponen- 
tiated give the finite phase-space translation operators, 
X(s) = e- zsp and Z(t) = e m , with s,t € R. X(s) acts 
on a continuously indexed computational basis state \q) q , 
an eigenstate of q, &sX(s)\q) q = \q + s) . Eigenstates of p 
transform similarly: Z(t)\p) = \p + t) . Transformation 
between the position and momentum basis is given by 
the Fourier transform operator F — exp[i(7r/4)(g 2 +p 2 )}, 
with F\s) = \s) p . This is the generalization of the 
Hadamard gate for qubits. The controlled operations 
CNOT and CPHASE are generalized to controlled- V 
(Cx) and controlled-Z (Cz), respectively. These opera- 
tors effect a phase-space displacement on the target by 
an amount determined by the position eigenvalue of the 
control: Cx = cxp(—iq®p) and Cz = exp(iq®q), where 
the order of the systems is (control (g> target). 

The essence of the qubit cluster-state model of QC lies 
in the one-qubit teleportation circuit 0, 0] . This cir- 
cuit gives the ability to teleport operations diagonal in 
the computational basis onto the state in question af- 
ter the cluster has been prepared. This allows dynamics 
to be performed solely through measurement. The CV 
analog of the one-qubit teleportation circuit is 
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In this diagram, |0) = (2n)~ 1 / 2 J dq\q) q is a zero- 
momentum eigenstate (the generalization of |+)), the 
controlled operation indicated is a Cz gate, and D is 
any operator diagonal in the computational basis (i.e., 
of the form exp[if(q)]). The projective measurement is 
of q and yields a real number s, which becomes the ar- 
gument of the displacement X(-) at the output of the 
circuit. The essential feature of this circuit is that the 
Cz gate commutes with any diagonal operator D. This 
means that even though D is applied after the Cz gate, it 
acts as if it had been applied before. Since the operations 
D and F^ followed by computational basis measurement 
are equivalent to a single measurement of D^pD, manip- 
ulating quantum information in the CV cluster is possi- 
ble through projective measurements alone. Concatena- 
tion of these circuits makes it possible to implement any 
single-mode unitary ||. 

As is the case for qubits [13], every CV cluster state 
has a graph state representation, where each node in the 
graph is a separate CV mode, and each link in the graph 
represents a Cz that has been performed between the cor- 
responding nodes (systems). Linear graphs, where each 
node has at most two links, can be used for single-mode 
evolutions, but not multi-mode gates. The simplest im- 
plementation of a Cz gate involves a graph state with a 



link between two adjacent quantum wires: 
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The lines to the left of nodes 1 and 2 indicate that a 
bipartite state \ip) will be teleported down two quantum 
wires to arrive at nodes 1 and 2. Measuring p on nodes 1 
and 2 leaves (X(st)F® X(s 2 )F)C z \tp) on nodes 3 and 4. 

A small set of Hamiltonians that are polynomials in q 
(e.g., {q,q 2 /2}), along with the Fourier transform, are 
sufficient to implement any single-mode Gaussian 
Furthermore, adding the ability to perform a Cz oper- 
ation (as described above) allows implementation of all 
multi-mode Gaussians. While this is not sufficient for 
universal QC, given an encoding that maps all qubit Clif- 
ford operations to CV Gaussian operations (the GKP en- 
coding being one example 9] ) , this would be sufficient for 
many quantum error correction protocols .21,]. Adding to 
the toolbox any single non-Gaussian projective measure- 
ment allows for universal QC using CV cluster states || . 

Optical Implementation — Since each mode of the 
electromagnetic field behaves as an independent har- 
monic oscillator, we can use these modes as CV systems 
for our CV cluster state. To do this, we choose the com- 
putational basis to be the "position" (amplitude) quadra- 
ture of quantum optics for each mode. The "momentum" 
(phase) quadrature for each mode becomes the conju- 
gate basis. The commutation relations [a, aS] = 1 and 
[q,p] = i are satisfied by the definitions q = (a + a))/\/2 
and p = —i(a — al)/>/2 for each mode. In this unitless 
convention, the variance of the vacuum state (which can 
be measured experimentally using homodyne detection) 
is given by (g 2 ) = (p 2 ) — 1/2. 

Construction of an ideal CV cluster state requires zero- 
momentum eigenstates, which cannot be normalized and 
are thus unphysical. In this optical model, they represent 
infinitely squeezed vacuum states, which require infinite 
energy. We can approximate them, though, by finitely 
squeezed vacuum states: 



\0,n) p :=(nfl 2 ) 
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with £1 2 < 1 being the variance of a Gaussian wave packet 
in momentum space (with (p 2 ) — il 2 /2). The states 
|0, f2) are defined analogously with p — > q in Eq. J2J. 
Note that |0,O) = |0,fi _1 ) . The fact that these states 
are finitely squeezed means that we will not have perfect 
fidelity while propagating quantum information through 
our cluster. This will be addressed later. Given the 
graph state that we wish to create, we need one indepen- 
dently squeezed mode per node, and we need the abil- 
ity to perform a Cz gate between modes in accordance 
with the graph. This operation is a quantum nondemo- 
lition (QND) interaction and can be implemented using 
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two beamsplitters and two in-line squeezers [22J. Alter- 
natively, it could be directly realized via an optical-fiber 
cross- Kerr interaction J2j| in the Gaussian regime of large 
photon number j^. (See also Sec. Ill of Ref. 01 f° r f ur ~ 
ther ideas.) 

Propagation down a quantum wire (D = I) is achieved 
through momentum-quadrature homodyne detection. As 
discussed previously, multi-mode Gaussian operations re- 
quire only that we can apply D = e ls ^ and D = e 4 ** I 2 
for all s,t£l. Applying a gate D to the encoded state 
is achieved by measuring the operator D^pD. Thus, the 
Z(s) — e ls ^ gate is implemented by measuring the op- 
erators Z(—s)pZ(s) = p — s. This is trivial to imple- 
ment: simply measure p and subtract s from the result. 
The gate denoted P(t) — exp(itq 2 /2) is implemented by 
measuring P(—t)pP(t) — p + tq. Notice, however, that 
by defining 9 = tan _1 (— t), we can rewrite this opera- 
tor as (pcosO — q sin 8)/ (cos 0), which is simply homo- 
dyne detection in a rotated quadrature basis, followed 
by a rescaling of the measurement results by a factor of 
cos0 = (l + < 2 ) -1 / 2 . Thus, once the cluster has been pre- 
pared, we are able to perform all multi-mode Gaussian 
operations simply through homodyne detection. 

Furthermore, analogously to implementing Clifford 
group operations on qubit cluster states, all multi-mode 
Gaussian operations may be implemented on CV clus- 
ters with the appropriate measurements made in any or- 
der. Performing the measurements in a different order 
is equivalent to commuting Gaussian operations through 
the (Gaussian) measurement-dependent corrections, re- 
sulting in different corrections, but leaving the measure- 
ment bases unchanged. This is known as parallelism in 
cluster-state QC [23. Non-Gaussian operations in gen- 
eral cannot be parallelized, since later measurement bases 
will depend on current measurement results, a property 
known as adaptiveness. 

Universal QC requires the ability to implement at least 
one non-Gaussian operation In our case, this will be 
achieved through a measurement in a non-Gaussian basis. 
While one can, in principle, use the continuous degree of 
freedom directly for QC, it will almost certainly be more 
practical (considering experimental errors) to encode fi- 
nite dimensional systems in the CV modes, e.g., as in the 
GKP proposal |S| , which encodes one qubit into each os- 
cillator. In this case, the optimal non-Gaussian operation 
would be tailored to implement a desirable non-Clifford 
unitary in the qubit space. Photon counting is one pos- 
sibility and fits nicely into the cluster formalism since it 
is already a projective measurement. Another option is 
to measure in a nonlinear polynomial basis, such as that 
corresponding to the observable p + uq 2 for any one par- 
ticular choice of u. This is equivalent, in the language of 
Circuit 0J, to implementing the gate D — e m * / 3 . The 
GKP proposal discusses both options in more detail. We 
leave the questions of encoding scheme and non-Gaussian 
measurement to future work. 



Experimental Errors — Possible sources of experi- 
mental error include the finite squeezing of the input 
states, mixed input states (but still Gaussian), and dis- 
tortions due to the QND operation used to form the clus- 
ter. Since any physical implementation of our protocol 
will be forced to use finitely squeezed states (because of 
finite energy requirements), we will consider the effects 
of finite squeezing in some detail. 

Finite squeezing in Eq. J2J modifies the output of the 
circuit in Circuit to M.X(s)FD\ip), where M. is a 
distortion that applies a Gaussian envelope in position 
space with zero mean and variance fl~ 2 : 

M\i>)0C J dqe^/ 2 \q) qq (q\^) . (4) 

Notice that this is not a unitary transformation, and the 
state must be renormalized after this envelope is applied. 
This is also equivalent to convolution in momentum space 
by a Gaussian with variance fl 2 . Mixed input states can 
be accommodated in this analysis (in the Wigner repre- 
sentation) simply by allowing the convolution width to 
be independent of the width of the Gaussian envelope. 
Thus, the transformation implemented by each measure- 
ment, which used to consist solely of Fs, Ds, and phase- 
space displacements, now includes a ubiquitous distor- 
tion at each step in the evolution. The severity of this 
distortion depends inversely on the amount by which the 
sources are squeezed. 

Concatenated circuits of the form Q apply the trans- 
formation ■ • • MX(s 2 )FD 2 MX(s 1 )FD 1 \ip) to the in- 
put. Alternatively, we can gather the fixed distortions 
to one end of the operation^and transform this into the 
useful form Uo(s\, ■ ■ ■ , s n )M(si, . . . , s„) \tp), where Uq is 
the unitary that would be applied in the case of an ideal 
cluster, and M.(s\, . . . , s n ) is a distorted initial state, 
with the distortion now depending on both the measure- 
ment results and the gate to be implemented. The ef- 
fect is the same for multi-mode gates: at each measure- 
ment step, a fixed distortion M. is applied to each mode. 
Specifically, in the case of the Cz gate, the resulting out- 
put is (MX( Sl )F <S> MX(s 2 )F)C z \f). The distortion 
operations in the multi-mode case can similarly be gath- 
ered to the right while becoming measurement- and gate- 
dependent. Errors in the QND operation can be modeled 
as additive Gaussian noise, which has a similar distorting 
effect, the strength of which scales as the number of links 
in the cluster's graph. 

Experimental Proposal for Cluster-Based Error 
Reduction — Parallelism, which is a feature particu- 
lar to cluster-state QC, along with post-selection, can 
be used to reduce the impact of the errors described in 
the last section when implementing Gaussian operations. 
We propose an experiment to demonstrate this. For con- 
creteness, consider a linear cluster of five nodes (although 
any number greater than three will suffice): 

-KTW2M3M3MI) (5) 
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(The line to the left of the first node indicates where this 
cluster might be attached to another one.) Many simple 
Gaussian operations may be implemented on this clus- 
ter through homodyne detection on the first four nodes. 
With each such measurement there is the possibility of 
the resulting distortion severely affecting the quantum 
state in a measurement-dependent way (see the previous 
section). However, if we choose to delay applying the 
QND operation that connects node 1 to node 2, we can 
isolate nodes 2-5 into a "mini-cluster," which is separate 
from the quantum state to be acted upon. By measur- 
ing nodes 3 and 4 before attaching the mini-cluster to 
the input state, we can calculate the effect of the dis- 
tortion from these two nodes before that distortion ever 
affects the state. If this distortion does not preserve the 
Wigner phase-space region likely to be occupied by the 
input state (which depends on the chosen encoding), we 
discard this mini-cluster and try again. If it does, we 
perform the QND operation to attach nodes I and 2. We 
now have only two "dangerous" measurements to make 
(on the newly attached nodes) instead of four, with the 
output appearing on node 5. State tomography can be 
used to compare the fidelity of these two approaches. 

This technique generalizes easily to multi-qubit oper- 
ations and can, in fact, be applied to mini-clusters im- 
plementing any Gaussian operation. The greatest ben- 
efit will be for Gaussians that require many measure- 
ments. While we have "bent the rules" of cluster-state 
QC a bit by delaying attachment of the mini-cluster and 
by post-selecting mini-clusters based on measurement re- 
sults, this may yet prove to be a practical procedure for 
dealing with experimental errors. This result has the 
flavor of Ref . |26| , wherein it is shown that through post- 
selection the reliability of an error-correcting ancilla clus- 
ter (called a "telecorrector" ) can be guaranteed before it 
is attached to the state to be corrected. 

Conclusion — We have generalized the notion of uni- 
versal cluster-state quantum computation to continuous- 
variable systems. We have proposed an optical imple- 
mentation that uses squeezed light sources and quan- 
tum nondcmolition operations to build a Gaussian clus- 
ter state. Homodyne detection alone suffices to imple- 
ment all multi-mode Gaussian operations using the clus- 
ter state, with the addition of one non-Gaussian measure- 
ment allowing for universal quantum computation. Many 
of the properties of qubit-cluster computation also apply 
to the continuous- variable case, including parallelism and 
adaptiveness. Within the continuous-variable approach, 
a lower-dimensional encoding scheme will most likely be 
required for experimental viability. Due to their Gaus- 
sian nature and deterministic method of construction, 
we expect that continuous- variable cluster states will al- 
low for further experimental demonstrations of the prin- 
ciples of cluster-state quantum computation. We have 
proposed such an experiment to demonstrate improve- 
ment in the fidelity of Gaussian operations using post- 



selection and parallelism. 
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